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Abstract 

Nuclear and lattice physicists share several topics of common interest, like hadronic 
masses, electroweak form factors and structure functions of hadrons. The main physics 
issues that a new collaboration between a nuclear physicist and the lattice group of 
the SPQ c dR collaboration is going to address, are summarized and preliminary results 
obtained for meson and baryon masses are presented. The importance of adopting non- 
perturbative approaches based on the fundamental theory of the strong interaction, like 
lattice formulations of QCD, is stressed. 
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1 From nuclei to quarks and gluons 



Since several years an increasing part of our nuclear physics community is deeply involved 
in the study of several topics pertaining to hadronic physics, like hadron mass spectra and 
decays, form factors of the nucleon and its resonances, structure functions of hadrons, quark- 
gluon plasma, etc. Since Quantum ChromoDynamics (QCD) is believed to be our funda- 
mental theory of the strong interaction, quarks and gluons are the degrees of freedom to be 
used in the description of the hadron structure. This fact implies two requirements. The 
first one is to describe the hadron structure through a non-perturbative approach based on 
the fundamental theory. The second one is that the same approach should be used to calcu- 
late different matrix elements relevant for QCD. Both requirements are fulfilled by lattice 
formulations of QCD. Thus it is clear that lattice QCD represents the future of theoretical 
nuclear physics. 

The aim of this contribution is to point out that there are several important physics 
issues shared by the nuclear and lattice communities. In this respect a new collaboration 
has recently started between myself and the lattice group of the Southampton-Paris-QCD- 
Rome (SPQ c dR) collaboration with the aim of addressing: i) the meson and baryon mass 
spectra; ii) the electroweak form factors and structure functions of the baryon octet and 
decuplet; iii) the electric dipole moment of the neutron induced by CP violations both in 
the electroweak and in the strong sectors (see in the latter case Ref. [[!]]). 

After recalling the main features of the Wilson and Clover formulations of QCD on 
the lattice our preliminary results obtained for meson and baryon masses in the quenched 
approximation will be presented. 

2 Lattice QCD 

Lattice formulations represent a possible way to regularize the ultraviolet divergencies of 
a quantum field theory. At the same time they are a very powerful tool to investigate 
nonlinear field theories in their own non-perturbative domain. As far as QCD is concerned, 
various lattice formulations have been proposed in the literature and are under continuous 
investigation and improvement. The SPQdR collaboration adopts the Wilson formulation 
of the QCD action in which the total action is given by the sum of a fermionic and a 
gluonic part 

Swilson = S$ + Sq , (1) 

which in Euclidean space read as 
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and 



S G = — \ 1 ~~ o R<3 t Tr Uplaquette] \ ■ (3) 

90 plaquette ^ ° ' 

where ipf is the fermionic field with flavor /, mj is the bare quark mass, go is the bare 
coupling constant, a is the lattice spacing and U^{x) = exp(ig aA c fM X c /2) is the gauge link. 
In Eq. (0) the terms proportional to the parameter r represent the so-called Wilson term, 
which avoids the fermion doubling problem. The Wilson term vanishes in the continuum 
limit (a — > 0), but breaks chiral symmetry. In Eq. (Q) the sum extends over the symmetric 
plaquettes, i.e. all the plaquettes having x as one corner with the plaquette operator given 
by the product of the four gauge links comprising the plaquette in the clockwise sense. 

In the continuum limit both and Sq reduce themselves to the fermionic and gluonic 
part of the QCD action, respectively. However, the convergence is different for and Sq, 
namely 

SV> -> [Sqcd)tp + °( a ) ' 

S G - [S QCD } G + 0(a 2 ) . (4) 

Therefore, a modification of the Wilson action has been proposed in order to kill all the 
linear terms in the lattice spacing. The new action is known as the Clover action || and 
reads explicitly as 

Sciover = Swilson ~ Csw<^ Yl '90J ^ f(x)a^ u P^ f(x) (5) 

x,f ^ 

where is the plaquette operator and c$w is a parameter which depends on g and should 
be determined non-perturbatively. This was done in Ref. 0] and one finally gets 

Sciover — * SqcD + 0(d 2 ) (6) 

Thus, the Clover action (|5]) makes the convergence to the continuum limit faster as 
fas as the action is concerned. However, the improvement provided by the Clover action 
does not exhaust the full improvement program, because in general fermionic operators have 
to be improved themselves in order to guarantee that their matrix elements are not affected 
by linear terms in the lattice spacing. Nevertheless, as far as hadronic masses are concerned, 
it is enough to consider only the improvement of the action. 

Hadronic masses can be determined by analyzing the time evolution of the two-point 
function at zero-momentum 

G H {t) = Y(0\T{O H (Z,t) O H (0,0)}|0) (7) 

X 

where T stands for the time-ordered product and On(x,t) is an interpolating field having 
the same quantum numbers of the hadron H. Inserting a complete set of hadronic states 
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one gets Guif) = J2 n Zn ex P{~Mnt), where is the mass of the n-th eigenstate with 
the quantum numbers of the hadron H and is the overlap of the interpolating field 
with the n-th eigenstate. Therefore, for large source times the ground state is isolated, 
i.e. Gn(t) —>t»a Zq exp(—M^t). The value Mq can be obtained from the plateau of the 
effective mass M e ff(t), defined as 

M eff (t/a) = \n[G(t - a)/G(t)} ^ t>>a M H • a (8) 



3 Meson masses 

We have carried out QCD simulations in the quenched approximation using 320 gauge 
configurations in a 24 3 x 64 lattice volume with a bare coupling constant go corresponding 
to (3 = 6/(7q = 6.0. We have considered various values of the hopping parameter K, which 
is connected to the (lattice) bare quark mass m by 



1 

m = — 
2a 



K K, 



(9) 



where K c is the critical value of the hopping parameter corresponding to the chiral point 
m = oi. In what follows SU (3)-flavor symmetry is assumed. 

In case of pseudoscalar (PS) and vector (V) mesons we have used the standard 
interpolating fields 

PS {x,t) = u(x)j 5 d(x) , 
Oy(x, t) = u(x) r y fl d(x) . (10) 

Taking into account the periodicity of the lattice the two-point function has the following 
large-time behavior [t » a and (T — t) » a] 

G{t)^Zo[e- Mot + e~ Mo( - T -V] (11) 

where in our case T = 64 • a. From Eq. (|TT|) one can see that the first and second halves of 
the lattice time contain the same information. The effective mass obtained for PS and V 
mesons at K = 0.13300 is reported in Fig. 1 to illustrate the quality of the plateau. 

The values obtained with the plateau method for the PS- and \^-meson masses are 
reported in Fig. 2, where it can be seen that our results are in very good agreement with 
those obtained in Refs. M, || at the same value of (3. The squared PS masses can be nicely 
fitted by a linear relation with the inverse hopping parameter [see Fig. 2(a)]. This means 
that Mp S is almost proportional to the lattice bare mass m. The intercept of the linear slope 



b The lattice mass m differs from the bare mass m° of Eq. (||) . Indeed, since the Wilson term breaks chiral 
symmetry, the quark mass receives an additive renormalization, i.e. in = m° — M , where M diverges in the 
continuum limit. Thus, using the Wilson (or Clover) formulation the chiral limit is not given by m° — > 0, 
but instead by m Q = M. 
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(b) V-meson J 
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Figure 1: Effective mass for PS (a) and V (b) mesons versus the source time t in units of 
the lattice spacing a for K = 0.13300 . The solid lines show the part of the plateau used to 
extract the ground-state masses. 



shown in Fig. 2(a) corresponds to K = K c = 0.13526 (2). From Fig. 2(b) it can be seen that 
the calculated V^-meson masses exhibit a linear relation with the squared PS'-meson mass, 
viz. 



M v -a = A + B (Mp S ■ a) 



(12) 



where A and B are two coefficients fitted to our lattice results. Then, using the physical 
masses of K and K* mesons we obtain for the inverse lattice spacing the value a~ l = 
1.90 ± 0.08 GeV, which corresponds to a = 0.104 ± 0.004 fm. Moreover, using the physical 
masses of the pion and the kaon the lattice bare masses can be determined from the slope 
shown in Fig. 2(a), obtaining mi = (m u + m^)/2 = 3.3 ± 0.1 MeV and m s = 78 ± 3 MeV. 
These masses have still to be renormalized in some scheme. Adopting the MS scheme one 
has m MS (n) = Z MS (afi) ■ m, where \x is the renormalization scale (conventionally taken at 
the value \i = 2 GeV) and Z MS is a renormalization coefficient which has been calculated 
nonperturbatively in Ref. g. At = 6.0 we get mf^(2 GeV) = 4.1 ± 0.1 MeV and 
m 1 f s {2 GeV) = 96 ± 4 MeV which compare favorably with world averages (cf. Ref. [§). 



5 




Figure 2: (a) Squared PS-meson mass versus the inverse hopping parameter; (b) ^/-meson 
mass versus the squared PS mass in units of the lattice spacing a . The solid lines are a 
linear fit to our results (squares). Dots and diamonds correspond to the results of Refs. |5|. [|, 
respectively. 



4 Baryon masses 



In case of baryons we have considered two interpolating fields, having the same quantum 
numbers of the proton and of its negative-parity partner (the Sn resonance), namely: 



O n +(x, t) 
O n -(x, t) 



V6 
1 

76 



£-abc 



u T a C^ b d b 



l5Uc 



(13) 



where C is the charge conjugation matrix. However, one has On+(x, t) = — 75O7V- (x, t)^, 
which implies a structure of the baryon two-point function different from the meson one [see 
Eq. (|TT])] . Indeed, for t » a and (T — t) » a one has 



G{ t) ^ l±Jl\z +e - M+t -Z^ M '^ 
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Thus, in the upper components the proton propagates forward while its negative partner 
goes backward. The situation is simply reversed in the lower components of the two-point 
function. Therefore, from the large-time behavior of the two-point function one can obtain 
simultaneously the ground-state masses of a hadron and of its negative-parity partner. Our 
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Figure 3: (a) Positive-parity nucleonic mass versus the squared PS mass; (b) ratio of negative 
to positive parity nucleonic masses in units of the lattice spacing a . The solid lines are a 
linear fit to our results (squares). In (a) Dots and diamonds correspond to the results of 
Refs. |5|, |(J , respectively. 



results are collected in Fig. 3 and compared (favorably) with those of Refs. [|, [| . Adopting a 
linear fit with the squared PS mass [as suggested by Fig. 3(a)], we obtain for the mass of the 
proton the value 0.98 ± 0.05 GeV and for the mass ratio M N -/M N + the value 1.65 ± 0.06. 
Both results compare very favorably with the experimental data M p = 0.938 GeV and 
M Su /M p = 1.63 from 0. 

In Ref. [[J the effects of S'[/(3)-flavor breaking on light-hadron masses have been 
investigated in details and it was shown that to a very good approximation the light-hadron 
masses depend linearly on the sum of the lattice bare quark masses. Therefore, using the 
slopes shown in Fig. 3 we can estimate the masses of strange hadrons, like the cascade S 
and the A baryon. Our results for various hadron masses are collected in Fig. 4. Note that: 
i) the interpolating fields (p~3| ) do not distinguish between A and £ baryons, and ii) the 
masses of tt, K, p and K* mesons are not reported, because they are used to fix the chiral 
point, the lattice spacing and the lattice bare quark masses. From Fig. 4 it can be seen that 
several experimental masses are nicely reproduced by our quenched calculations with the 
only remarkable exception of the A (1405) resonance. The anomalous light mass of the latter 
suggest a possible interpretation as a NK molecule, which is totally missed in a quenched 
calculation. 

The results presented have been obtained at (3 = 6.0 only. Therefore, QCD simu- 
lations at larger values of [3 (which correspond to lower values of the lattice spacing) are 
needed to assess the uncertainties related to the continuum extrapolation. We stress that 
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Figure 4: Mass of various light hadrons. Open squares are our lattice results, while full dots 
represent experimental data from PDG 0. 



thanks to the use of the Clover action (|5|) the above-mentioned uncertainties are only of the 
order of 0(a 2 ). 



5 Conclusions 

Several important physics issues are shared by nuclear and lattice communities. In this 
respect a new collaboration between a nuclear physicist and the lattice group of the SPQdR 
collaboration has recently started with the aim of addressing: i) the meson and baryon 
mass spectra; ii) the electroweak form factors and structure functions of the baryon octet 
and decuplet; iii) the electric dipole moment of the neutron induced by CP violations both 
in the electroweak and in the strong sectors (see for the latter case Ref. |IJ). Preliminary 
results obtained for meson and baryon masses in the quenched approximation have been 
presented. 

Our main conclusion is that the future of theoretical nuclear physics is represented by 
lattice QCD, to which efforts and manpower should be devoted by our nuclear community. 
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